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Abstract 

We study the problem how to deal with tensor-type two-loop integrals in the Loop Regularization 
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find a new divergence structure: the regulated result for each two-loop diagram contains a gauge- 
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divergences cancel, recovering the gauge invariance and locality. 
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I. INTRODUCTION 



This is the second of a series of papers to study the consistency of Loop Regularization 
(LORE) scheme and the general structure of QFT models in the higher-loop calculations. In 
our first paper 01, we have shown that the LORE method is consistent for general scalar-type 
two-loop integrals. We studied the general q;/37 integralsjsf and mainly solved the problem of 
disentangling the overlapping divergences, which is the main new feature of the perturbative 
calculations beyond one-loop order. We argued that the overlapping divergence structure 
would, in general, transforms some of the ultraviolet (UV) divergences into the integrals 
of the ultraviolet divergence-preserving(UVDP) parameters (or Feynman parameters). We 
found that the LORE method was particularly suitable to regularize the UV divergences of 
this kind since we can easily transform the UVDP parameter integrals into the irreducible- 
loop-integral(ILI)-like integrals by multiplying the integration variable by a mass scale, which 
is exactly the object regulated in the LORE method. Furthermore, it was very useful 
to introduce the Bjorken and Drell's analogy between the Feynman diagrams and electric 
circuits, which allowed us to identify the UV divergences of the UVDP parameters with 
those contained in the subdiagrams of the original Feynman diagrams. Then we applied 
all the tools we found to the study of two-loop calculations of the 0^ model and obtained 
the consistent results for /3-function of the coupling constant and mass parameter. Finally, 
we gave the general procedure of the application of the LORE method to the higher-loop 
calculations, even beyond two-loop order. 

In the present paper, we will focus on one of the important remaining problems- how 
to deal with the tensor-type integrals at higher-loop order. Recall that at one-loop order 
in LORE tensor-type ILIs and scalar-type ILIs can be related by introducing consistency 
conditions, which relate these type of integrals with the same divergent behavior. So the 
immediate question is whether these consistent conditions, which should be applied to the 
1-folded ILIs of each loop momentum variable, is still consistent with the general features we 
find in our previous paper. In particular, we want to ask whether the consistent conditions 
can guarantee the cancellation of various harmful divergences. Another question involves 
the gauge invariance which is ensured by the consistency conditions in LORE as shown in 
0, Hf. Thus, it is necessary to inquire whether the consistency conditions can still preserve 
gauge invariance and the associated Ward identities. In order to study these questions, we 
use the massless Quantum Electrodynamics (QED) as our simplest example, which involves 
the tensor-type integrals and the gauge invariance at the same time. Specifically, we shall 
calculate the two-loop photon vacuum polarization diagrams in detail. As will be shown in 
the following calculation, the LORE method, with the prescription of consistency conditions, 
can give the sensible results which satisfy all the requirements: gauge invariance and locality. 
We should emphasize that the result above only requires the use of consistency conditions 
found at one-loop order and does not need to introduce additional two-loop order consistency 
conditions. 

One special feature worth mentioning is the appearance of a new divergence structure: 
each two-loop Feynman diagram contains a quadratically harmful divergence even when we 



combine the result with the corresponding one-loop counterterm insertion diagrams. These 
new divergence structure is expected and compatible with the usual power counting for each 
diagrams. However, these quadratically harmful divergences will cancel with each other when 
we sum up all three group of diagrams in this order, rending the final result local. Another 
way to see this cancellation can be obtained from the requirement of gauge invariance of 
the final result since these quadratically harmful divergences break the gauge invariance and 
must disappear in the final result. We argue that this is a general feature for gauge theories, 
since quadratic divergences would in general give the photon a mass, which is incompatible 
with gauge invariance. The appearance of this new divergence structure explicitly shows one 
of the main advantage of LORE method that the LORE method enables us to calculate the 
full result for any Feynman diagram, especially the quadratic divergences, which is absent 
in many other regularization methods, like the popular dimensional regularization (DR). In 
the light of detailed and complete calculations, we conclude that the LORE method can 
properly regularize tensor-type high-loop integrals (at least up to two-loop) and consistently 
apply to general gauge theories, like the Standard Model (SM) of particle physics. 

The paper is organized as follows. Section 2 gives a short introduction to the LORE 
method. In Section 3, we first compute all of the one-loop one-particle- irreducible diagrams 
that need to be regulated, in order to set our notations. Then we present the results of the 
massless QED vacuum polarization at two-loop order by using the LORE method, leaving 
the many of the details to the Appendix B and C. Appendix A contains some useful formulae 
in the LORE method. Appendix D contributes to the detailed calculation of a new doubly- 
logarithmic divergent integral encountered in our two-loop calculations. 



II. IRREDUCIBLE LOOP INTEGRALS(ILIS) AND THE PRESCRIPTION OF 
LOOP REGULARIZATION (LORE) 

In this section, we give a brief introduction to the Loop Regularization (LORE). A more 
detailed discussion will be found in the seminal paper by one of the authors 

It has been shown in0, [sf that all one loop Feynman integrals can be reduced into the 
following 1-fold ILIs by using the standard Feynman parameterization method: 
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with I2 and Jq corresponding to the quadratic and logarithmic divergent integrals. Here the 
effective mass factor M"^ is a function of the external momenta Pi, the masses of particles 
rrii and the Feynman parameters. 

When the regularized 1-fold ILIs satisfy the following consistency conditions 0, [sj: 

tR tR tR / \ tR 
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the resulting loop corrections are gauge invariant. Here the superscript "R" denotes the 
regularized ILIs. 



Note that the introduction of the irreducible loop integrals (ILIs) is crucial in the loop 
regularizationj^, [El, and it has been shown that all Feynman loop integrals can be expressed 
in terms of the ILIs. In the definition of ILIs, one important feature is that there should be 
no factor of k"^ in the numerator of loop integration and all the ILIs can be classified into 
the scalar type ILIs with the following loop integrand 
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and the tensor type ILIs with the following loop integration 

kfxhjy * • • hp 

In manipulating the Feynman loop integrals into ILIs, one should always perform the 
Dirac algebra and Lorentz index-contraction first to obtain the ILIs defined by the above 
"simplest" forms for the one loop case. For example, the integrand 

■ k,K/{e - M^f 

should not be written as 
Instead it should be expressed as 

- My 

and then rewrite the k"^ in the numerator in the form (A;^ — M^) + so as to cancel out 
the first term by the denominator to get, 

g^" ■ k,k,/ie - = h + M^- Jo. 

A simple regularization prescription for the ILIs was realized to yield the above con- 
sistency conditions. The procedure is: rotating to the four dimensional Euclidean space 
of momentum, replacing the loop integrating variable k"^ and the loop integrating measure 
J d^k in the ILIs by the corresponding regularized ones [k'^]i and j[d'^k]i: 



e ^ [k% = k^ + Mi , 
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J d^k^ J [d% = hm J^cf J d^k (3) 

where Mf (/ = 0,1, ■ ■ ■ ) may be regarded as the regulator masses for the ILIs. The 
regularized ILIs in the Euclidean space-time are then given by: 



tR 

-'-2a 



N,M^ j-^ J 



" d^k 1 



' 1=0 

N 

R 



(27r)4(A;2 + M2 + M2)2+°' 
r (j^u k k k k 

^-2a^.upa H ^™2Z^^' J (27r)4 + M2 + M2)4+" 



where the coefficients cy are chosen to satisfy the following conditions: 
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with the notation lim^y^^fa denoting the limit limjv,M2-5>oo- may take the initial conditions 
Mq = /i^ = and = 1 to recover the original integrals in the limit oo (/ = 1, 2, • • • 

). Such a new regularization is called as Loop Regularization (LORE) [3,[5|. The prescription 
in LORE method is very similar to Pauli-Villars prescription, but two concepts are totally 
different as the prescription in the loop regularization is acting on the ILIs rather than on 
the propagators as in Pauli-Villars scheme. This is why the Pauli-Villars regularization 
violates non-Abelian gauge symmetry, while LORE method can preserve non-Abelian gauge 
symmetry. 

A simple solution of eq. (jSD, is to take the string- mode regulators, 

= /x^ + iMl (6) 
with / = 1, 2, ■ ■ ■ , and the coefficients cf to be of the form. 



Here Mr may be regarded as a basic mass scale of loop regulator . It has been shown in [5| 
that the above regularization prescription can be understood in terms of Schwinger proper 
time formulation with an appropriate regulating distribution function. 

With the string-mode regulators for and cf^ in above equations, the regularized ILIs 
I2 and Iq can be calculated to give[3, 

= T=^AMl - /.^[In ^ - 7. + 1 + 
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with fi^ = + M^, and 
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This indicates that the fig sets an IR 'cutoff' at = and M^ provides an UV 'cutoff'. 
For renormalizable quantum field theories, M^ can be taken to infinity (Mc — )■ 00). In a 
theory without infrared divergence, fig can safely run to /i^ = 0. Actually, in the case that 
Mc — )■ 00 and fig = 0, one recovers the initial integral. Also once M^ and are taken to 
be infinity, the regularized theory becomes independent of the regularization prescription. 
Note that to evaluate the ILIs, products of 7 matrices involving loop momentum such as 
//7^^ should be reduced to one of the independent components: 7^, 0"^,^, 757^^, 75 without k. 



III. MASSLESS QED VACUUM POLARIZATION AT TWO LOOP ORDER 



In this section, we will present our results of vacuum polarization in the massless Quantum 
Electrodynamics (QED) at two-loop order. The motivation for this computation is two- 
folded: (1) to show how to apply Loop Regularization (LORE) method to tensor-type two- 
loop integrals; (2) to give the first explicit example of applying LORE to the realistic model. 
We will show that LORE can preserve Ward Identity and in turn the gauge invariance of 
QED at two-loop order. Also, our result reproduces the results in the standard textbooks like 
[6-8]. It is an essential step towards applying LORE to the general gauge theory calculations, 
such as the Standard Model of particle physics. 

The Lagrangian of the massless QED is: 

CQED = -\{F^.f + i^i:^^ (10) 

where, = <9^ + ieA^ is the covariant derivative and F^^ = d^A^, — d^Afj, is the field 
strength tensor of electromagnetic^ field. The detailed Feynman rules for QED are referred 
to the standard texts, such as 



A. Regularization and Renormalization of QED at One-Loop Level 

In this section, we present the regularization and renormalization for all the divergent one- 
particle-irreducible (IPI) diagrams in the massless QED at one-loop level with the LORE 
for use in the later calculation. At one-loop level, the massless QED has three divergent IPI 
diagrams requiring regularization and renormalization. They are shown in Fig.(IT]) 
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FIG. 1: 1-Loop IPX diagrams: Left: Vertex correction; Middle: Electron self energy; Right: Vac- 
uum Polarization 

Using the standard method of Feynman parameters, we get for the vertex correction 

^ 2r(3)e3 1^ - - -)'^- " " ' 

where M^^ = {xpi — yp2Y — xpf — yp\. We can write this in terms of 1-folded ILIs as, 
_,rM(i) = 2r(3)e=^y" dxdy{{r^^'^^)llp-[{l-x)Tj,, + yj^2mxj^i + {l-ymi-2} 

= 2r(3)e^ j dxdyi—^^-I^ - [(1 - x)^i + y^2mxjfi + (1 - ymi-2. } (12) 



where in the second hne we have apphed the consistency condition Eq.([2]) to deal with the 
regularized logarithmic divergent integral 1^^^^. By using the explicit expression of 1^ Eq.® 
and integrating out J_2, we can write our result of the one- loop vertex correction as, 

Stt^ J /if [xpi - yp2) - xpi - 

where /xf = fi1+M? p^. According to our renormalization scheme in LORE method discussed 
in our first paper j^], we should define the one-loop vertex counterterm as: 

- le-ff'Si = -^1^ / dxdyihi ^ - 7^ = - In ^ - 7^ , 

where /i is introduced as the renormalization scale. With such definition of the renor- 
malized vertex correction at one-loop level is: 



8vr^ J ni {xpi - yp2y - xpi - ypi 
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where we have taken the limit — > 00 so that yo{-^) — ^ to simplify our final result. 

For the electron self-energy diagram, we can use the same procedure while the calculation 
is more straightforward. 
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Using the general formula of LORE, we can regularize the above logarithmic divergent 
integral as follows: 
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If we introduce the counterterm for the self-energy as: 

= -— ^^(In^ 
then the renormalized self-energy for electron is: 



mn^P) = / dx{l - x) In ^ — , (15) 

Note that the counterterm coefficients defined above have the interesting relationship = (52, 
which is the result of Ward identity in QED. This relationship is crucial to guarantee the 
consistency of our later discussion of renormalization at two-loop level. 



Finally, the one-loop vacuum polarization diagram can be written as: 
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By computing the trace of gamma matrices, we obtain the following expression in terms of 
ILIs: 

= -4e2 / (ix{(2/r'^ - g^"!^) - 2x{l - x){p^p'' - g^'p")!^]. 
Jo 

With help of consistency condition Eq.([2]), it is seen that the quadratically divergent terms 
in the first parenthesis cancel each other exactly, and the remaining term is proportional to 
{p'^p'^ — g^^p^), as required by gauge invariance. Using the formula for Iq, we get. 

Again, if we choose the counterterm for the photon vacuum polarization as: 

iip^'p" - g^^P^s = TT^ip^p'' - g^'p') / rfxx(l-x)(ln^-7j 

(pV-^7'^V)(ln^-7.) (17) 



the renormalized photon vacuum polarization at one-loop level is: 

-M^'^ = ^(P^'^ - ^^V) rfa;x(l - x) In ^ — , 

/7r"= Jq jjLg — x[l — xjp"^ 

where we have taken yo{ '^''~^^2^^^ ) when — )• oo as usual. This completes our 
discussion of the regularization and renormalization of massless QED at one-loop level. 



B. Self-Energy Insertion Diagrams 

Now we compute the photon vacuum polarization at two-loop order. In this subsection, 
we calculate the diagrams (ai) and (02). It is easy to see that the two diagrams are equal 
and can be recognized as the insertion of one-loop electron self-energy into the one-loop 
photon vacuum polarization. 

We first write 1-loop electron self-energy using the UVDP parametrization, rather than 
the Feynman parametrization like we did in Sec JIII A] The reason lies in the fact that this 
time the one- loop electron self-energy appears as the subdiagram in (ai) and may contribute 
to the overlapping divergence in the following calculation, 

2 duu r d'^k2 Iji 
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FIG. 2: Self-Energy Insertion Diagrams. Left: (ai); Right: (02) 

By inserting the above result into the expression of one-loop photon vacuum polarization 
diagram, we obtain: 



4 r du(u + 1) f d^h f d^k2 f\ 
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Note that we have used the usual Feynman parameter x to combine the denominators 
involving ki only, rather than combine all of the three factors once at all. This transforms 
the Feynman integral A^('*i) into the a/37— like integrals (we have already combine two 
factors in Eq. (|T8l) with UVDP parameter u ). The advantage of this procedure is that it 
effectively separates the UV divergences from the IR ones in the parameter space. To see this 
more explicitly, our first paper [Ij already explicitly showed that the a/37 integrals contain 
the most general UV divergence structure. In other words, the UV divergence can only be 
contained in a/37 integrals and their corresponding UVDP parameter integrals, rather in the 
Feynman parameter space. If we find some divergence in integrating Feynman parameter x, 
then it must be the IR divergence, rather than the UV one. Since we have already written 
the integral into a generalized a/37 form, for the rest factors, we will apply the UVDP 
parameters to combine them. 
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Note also that, in the above calculation, we used the traditional technique of completing 
squares for each internal loop momentum and then shifting the origin of it in order to 
eliminate the terms linear in the momenta. The resulting formula appears less symmetric in 
the UVDP variables than that with the general formulae in Ref . [2| . However, this traditional 
method is more flexible and more natural in practice when the numerator is a complicated 
function of internal and external momenta like the present case, rather than a constant for 
the scalar integrals. In the following calculation of Fig. (jl]), we will also use this traditional 
method. 



Since there is no cross terms between ki and k2 in the denominator, we can easily integrate 
out ki and with the LORE method sequentially: 
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where ij? ^ fi^— ^"^^^ (1 — ■f^)p'^ and A^2"o) denotes quadratic (logarithmic) divergent part 
from the integration of fc2- 

Note that since our integrals in the present paper do not contain any IR divergences, we 
can set /^^ — )■ which also plays the role of IR regulator in the LORE method. Moreover, 

2 2 

because the yi{-^) — )■ in the limit of — )■ 0, so these functions do not contribute to our 
final results and can be ignored. 

To make our discussion more concise, we only present our final regulated result for the 
diagram (ai), while relegating our calculational details of the UVDP parameter integrations 
to the Appendix iBl 
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(22) 



In order to see the cancellation of harmful or nonlocal terms for diagram (ai), we need 
to compute its corresponding counterterm diagram (a'j^), as shown on the left in Fig. ()3]), 
and this diagram contributes 





FIG. 3: Counterterm Insertion Diagrams for Diagram Group (a). Left: (a'^); Right: (og) 



Here we recognized that the internal momentum integration is essentially the same as the 
photon vacuum polarization at one- loop order and the result Eq. (fT7|) can be directly applied. 
By summing up the diagrams (oi) and {a[), we can obtain: 

[IGtt'^)^ —q^ 6 —q^ 36 

+ (^?^V - P'pni-k^r, ^ - 7.)' + ^(In ^ - 7.) + ^ In ^ In 
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1,2/^ 5 , /i^ 11 , yu' 1 17 TT^ 

— In — In — - H In — 1 \r. \^^) 

6 -ql 18 -p2 9 -g2 5 12 ig^^ ^ ^ 

Note that there is still a harmful divergence ^-^^t^^ g^'^M^ ln(— g^), which also breaks the 
underlying gauge invariance and its associated Ward identity. From the general argument 
of gauge invariance, we expect that terms like this should not appear in the final result of 
the two-loop vacuum polarization. As we will see in the following, the result of diagram 
(b) contains same term with the opposite sign. Therefore, the addition of (ai), (02) and 
(6) will eliminate this unwanted term, and recover the gauge invariance and locality of the 
underlying field theory. 



C. Vertex Correction Insertion Diagram 

Now we compute the more challenging diagram in Fig 4. Following the standard Feynman 
rules of QED, we can write down the explicit expression Ai^''^ 




FIG. 4: Vertex Correction Insertion Diagram (6) 

= (_,e)4 fd^f ^^_^)tr(Y-Y ' Y - 7- ' 

^ ^ '"W (2vrry (2vr)^^ '^'"^^ H. + H,^ Hi + H2+i Hi+f ' 

. 4 f d% f d% [\ r\ 

= te - — -r / - — — / axi / ax 



'-2 

tr[Yih - xii)Y{Hi + H2- xii)Y{Hi + ^2 + (1 - xi)^)7p(^i + (1 - 







{kl + Xi{l-Xi)p^f[{h + k2 + {x2-Xi)py+X2{l-X2)p^Ykl ^ ' 

where we have used Feynman parameters to combine some of the denominators. Next, in 
order to integrate over momentum ki, we apply the UVDP parametrization to combine the 
factors involving ki in the denominator to get 

1 

{kf + xi(l - xi)p'^y[{ki + k2 + {x2 - xi)py + X2{1 - X2)p'^Ykl 

- V(A\ A ~ ^'=1 ij^^ (^p^ 



where we define 

l + Vi 1 + f 2 

= (^1 + T^k, + "^^pf 

1 + f 2 1+V2 
(1 + fl)(l + 1^2) 1 + fl l + f2 

Then we make the following translation 

k,-^k,- -^^2 - ^^^P. (28) 

1 + t;2 1 + t;2 

By expanding the trace of the product of gamma matrices and contracting the Lorentz 
indices, we write the numerator as, 

= i-8){9^''k^i + [^k^ik^h ■ k - g^^klh ■ k - '^k^^llh ■ ki - 2l^k'(h ■ ki 
+2(7^"/3 ■k,h-ki + I'il^kl + I'^^llkl] + [I'ir^kl - It^llkl + g^^'h ■ kkl] 
+ [-2k';k'[k ■ h + g^'^klh ■ h + 2A;f/^/2 ■ ki + 2fcfZ^/3 ■ k^ - l^F^kl - l^F^kl] 
+ [-l^likl + Itr^kl + g^'^h ■ kkl] + [-2k>tklh ■ k + g^^klh ■ k + 2l^^klh ■ ki 
+2l'ik'(k ■ ki - rjlkl - l^l^kl] + [2k'^k'[k ■ k - g^^klh ■ k - 2l(^k'(h ■ ki 
-2k^^l\k ■ ki + 2g''''h ■kik-ki + l^^H^kf + l^l^^kl]} + 1'^!^^^, (29) 

where 1^1234 = t^^il^hl^hl'^hlph]- Note that in above expression, we have used ZjS to 
simplify our notation which represent the four factors in the Eg. (1251) after the translation of 
Eq. fl^Hl) . Due to their complexity, we do not write them out explicitly. 

After the manipulation above, the integration over the loop momentum ki is straightfor- 
ward, 

^ + :n^ii'j2 - + g'^'k ■ k) + TT^ig'^'k ■ k - m - m 



- ^3^1 + a'-'k ■ k) + TTv^ig'-'k ■ k - m - m 

Now we integrate over the loop momentum ^2- By using the UVDP parametrization, we 



combine the two factors containing A;| in the denominators 
1 du 1 



oo 



du- 



1 1 

^^r,, I 1 12 \(U I \2 212' (^-'-) 

~ , 1 1 

du- 



[rf'^fci], r(3) 



dU; 



(2vr)4 r(2)r(l) io [A:?-«A;i-M| 13 



k2,pi 

' 1 1 

''U (l+fi)(l+'y2) ^ l+«(l+i;i)(l+i;2)-^^ 

where we have defined 

2 _ r 1 r ^ Xi(l-Xi) ^ X2(l-X2) 

''o, J 

(l+Kl) (1+1)2) 



1 (X2 - Xi) 



2 



After performing Wick rotation of k2, we can integrate out k2 with the LORE method, 
and this naturally separates the Ai^^^ into four parts according to the different powers of 
the factor ^ i . The resulting expression is lengthy, so in order to keep the paper 

"+(l+«l)(l+i,2) 

in a readable length, we do not write it here. Appendix O gives the explicit form for each 
part, as well as the careful calculation of the remaining UVDP parameter integrations. We 
recommend the interested reader to resort to the Appendix O for such details. Below we 
only present our final results for the vertex correction insertion diagram (6): 

M^'^ ~ -7^2{9''M!{\n ^ - 7.) + ^^?^y (In ^ - 7.) 

Hg'^p' - P'PI ■ [^(In ^ - 7.) (In ^ - 7.) + ^(In ^ - 7.)] 
3 —p^ — 18 —p'' 
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-3 (In - 7.) + 3 (In ^^^^^ - 7.) " 2 " + r^^^' 

The counterterm diagrams {b[) and (62) for (6) can be computed directly and the result 

is 
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(34) 



FIG. 5: Counterterm Insertion Diagrams for Diagram (6). Left: {b[); Right: (62) 



and the summation of (6), {b[) and (62) gives: 

1 A/f2 fi A/f2 1 ,,2 ,,2 

Hg'V - P'Pn ■ [-^(In ^ - 7.)^ + ^(In ^ - 7.) + ^ In ^ In . ^ 



fi'^ 9 /i^ 3 — p2 — 

1 2 /i' 19 7 1 

6 —q^ 18 — 18 —q^ 6 

As expected from the argument in the end of the subsection IIII Bt the results of summing 
over the diagrams {b + b'l + 62) also possess a harmful divergence which violates the gauge- 
invariance . As we will see, this term has the right value to cancel the harmful divergence 
in the diagram group (a). 



D. The Results 



With the results Eqs.( l24l) and ( 135|) for the diagrams (ai + a'i) and (6 + 6']^ + 62) at hands, it 
is straightforward to obtain the final result of QED vacuum polarization at two-loop order: 

^final ^ 2A^(''i"'""'i) + 

:{9''Ml{\n ^ - 7.) + g'W^i^^ ^ - 7.) " I] 



{IQir'^y —ql 6 —ql 36 

Hrp' - P'f) [-^(In ^ - l.f + T^(ln ^ - 7.) + ^ In ^ In ^ 
6 /i^ 18 LL^ 3 —V —"^ 



/i^ is yU^ 6 —p^ —q^ 

1, 2 5 , /i^ 11 , /i^ 1 17 TT^, 

- In — ^ In — :r H In — \ \\ 

6 -g2 18 -p2 ^9 _g2 g ^ 12 18^^ 

8ie^ r „i.,.2/i Ml , 1 2/1 Ml 

{g^'MlilYi ^ - 7,) + -(7^V(ln ^ - 7^) 



(167r2)2 -ql '-^ ' g^ ^ ^ _g2 

+ (^'^V - pV) ■ [-^(In ^ - l.f + ^(In ^ - 7.) + ^ In ^ In 



jJ? Q jJ? 3 — p2 — g2 

1,2/"^ 19 , /i^ 7 , fi^ In 

- In^ H In H In aj] 

6 -g2 18 -p2 18 -g2 g ^JJ 



.(^My _ (In ^ _ + (In ^ - 7^)]. (36) 



(167r2)2 3 — p2 g _g2 

Note that there is a free mass scale —ql in the above expression, which was introduced 
to transform the divergent UVDP-parameter integrals into ILI-like one and the latter is the 



right object to be regularized in the framework of LORE method. If this scale is left in the 
final result, we must ask what is the physical meaning of this new scale and what the exact 
value is. In fact, as argued in our previous paper [1], this scale should be determined in terms 
of the intrinsic energy scales in the theory, like the external momenta or particle masses. 
In the present case, the only energy scale is the square of the external photon's momentum 
p^, which, in the massless electron case, needs to be space-like. Thus, if the —q^ has some 
physical meaning, it must have some relation to p^. The simplest possibility is that they are 
the same, —q^ = —p^. By applying the relation —q^ = —p^, we can reach our final result of 
photon vacuum polarization diagrams at two loop order: 

Q„-„4 -| T\/r2 

= 4;^(P>^-^^V)(lnz^-7.), (37) 

2 

where a = is the fine structure constant. 

In the light of our result, we have a few comments: 

1. As we have already seen, the quadratic divergences are canceled, and this is crucial to 
guarantee that photon does not obtain a mass from quantum fluctuation and that the 
whole theory remains gauge invariance. 

2. The harmful divergences like logp^ ■ logM^ and ■ logp^ vanishes also, which is 
expected as these terms are nonlocal and cannot be eliminated by any counterterms in 
the original Lagrangian which are local. Especially, the quadratic harmful divergence 
shows a new divergent structure, which can be expected from simple power counting 
in each of these diagrams. The nontrivial thing is that for every diagram this term will 
persist even when we include the corresponding counterterm insertion diagrams which 
only cancel the logarithmically harmful divergence. Only when we sum over all the 
relevant diagrams do the quadratic harmful divergences cancel each other, rendering 
the final result local. 

Moreover, as emphasized in the previous subsections, this feature is also expected from 
gauge invariance and the Ward identity in QED. Note that this quadratically harmful 
divergence is only proportional to g^^ and can give the photon a mass term, both of 
which are forbidden by the gauge invariance. We think this is a manifestation of the 
general phenomenon in gauge theories like QCD that although each diagram contains 
the gauge- violating quadratic divergences, but in the sum of all the relevant diagrams 
they cancel. The only difference is that this phenomenon occurs at two-loop order for 
QED while for QCD it happens at one- loop level already 0, Hj. It seems that it is the 
gauge invariance that rescues us from the harmful divergence "disaster" . 

3. We can compare our calculation with the LORE method with those using the dimen- 
sional regularization (DR). In the calculation with DR we cannot see such cancellation 
of quadratic harmful divergences because in DR, the final result comes from expansion 
around the actual space-time dimension 4 and in such a expansion only the logarith- 
mic divergences can be preserved while the quadratic or more higher divergences are 
simply omitted. It is the resulting incomplete expression in DR that makes this new 
divergence structure and its elimination invisible. However, since LORE enables us to 
calculate the full results in the exactly four dimensions, we can easily overcome this 
problem faced by DR, allowing the appearance of this nontrivial divergence structure. 



4. The final result Eq. (l37|) is transverse, which implies that our calculation respects gauge 
invariance and the corresponding Ward Identity. Therefore, the final divergence can 
be canceled by a single two-loop local counterterm for photon vacuum polarization. 

5. Our result agrees with other authors' results 



Therefore we come to our conclusion: LORE method, together with the prescription of 
consistency conditions, can be safely applied to general gauge theories, like QED, at least, 
up to two-loop order. 



E. Counterterm Diagrams 

Before ending, we want to show an interesting relationship between one-loop countert- 
erm insertion diagrams for photon vacuum polarization in massless QED. This can also be 
regarded as the check of the consistency for our calculations. We will prove that the coun- 
terterm insertion diagrams {a'l), (02), and (62) are, in fact, cancel each other due to the 
Ward Identity at one-loop order, though we have already given the results of these diagrams 
to show the cancellation of gauge-invariant harmful divergences in the proceeding subsec- 
tions. Obviously, A^("'i) = ^^("2) and Ai^'^'^^ = A^^'^a) since there is the up-down symmetry 
between {a[) and (02) as well as the left-right symmetry between {b[) and (62). Furthermore, 
we can write down the expression for Ai^""^^ and Ai^''^^ explicitly according to the ordinary 
Feynman rules: 

= / (38) 

(39) 

Adding up these four diagrams gives: 

(40) 

According to our renormalization scheme for the coupling constant and the electron wave- 
function at one-loop order, we have: 

5, = 62. (41) 

With the above result, the total contribution of the counterterm diagrams to the photon 
vacuum polarization vanishes. 



IV. CONCLUSION 



In this paper, we mainly discuss the problem how to deal with the general tensor-type 
high-loop integrals in the framework of LORE. We show that LORE, together with its 
consistency conditions for 1-folded ILIs, are enough to regularize the tensor-type integrals 
properly and can be consistently applied to the general gauge theories. In order to show this 
conclusion, we use the two-loop massless QED vacuum polarization as the simplest gauge 
theory example to show the general procedure. The final result is sensible: the two-loop 
correction is transverse and the harmful divergences and quadratic divergences are naturally 
canceled, which are required by the locality and the gauge invariance and its associated 
Ward identity of the underlying QED. Consequently, our result agrees with other authors' 
well-established one. Therefore, the computation of massless QED vacuum polarization at 
two-loop order explicitly shows LORE is consistent with the gauge invariance and give us 
confidence to apply it to other realistic theories, like the standard model in particle physics. 

What we want to emphasize is that from our present calculation, we have shown a new 
divergence structure: either the self-energy insertion diagram (ai) and (02) or the vertex 
correction insertion diagram (6) may contain gauge non-invariant quadratically harmful di- 
vergent term even when combined with their own one-loop counterterm diagrams. However, 
only when we sum over all diagrams in this order the harmful divergences cancel and the 
final result recovers the gauge invariance and locality. We argue that this feature is general 
for gauge theories, because the quadratic divergences which violate the gauge invariance 
will give the gauge boson mass and should be absent in the physical gauge invariant results. 
Furthermore, the reason why LORE method, rather than other regularization methods like 
dimensional regularization, can show this new divergence structure lies in the fact that 
LORE enables us to calculate the complete expression for any Feynman diagram, especially 
including the quadratic divergences, which is one of main advantages of LORE. 
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Appendix A: Useful Formula for the Regularization ILIs in LORE 

In the following, we will list some of the most useful regularized ILIs for reference: 



r 74 7 1 1 2^(0; — 2) 1 / A*^ X 1 

(A4) 

/ 1"*'*^!' (FWTW = - 1" f - 7» + 1 + 1} (A5) 

r 4 , kuK 1 . 2 ^(^-2) 1 ^ -I 

y ^]'(F + Mf + A1^)«+i ^ 2^''^'' r(a + l) (;.2).-2{ 1 " ?/-2(a-2)(^) } 

(A8) 

/ ^'"'^^WTW^' ^ ^ W}-^{^c^ - /.^[ In ^ - 7. + 1 + ]} 

(A9) 

/ [-^^^]- (.^+M^';';i^)^ - i w>^^^{ 1 - .-.(^) } (All) 

/r 74 7 ] kfikukpk„ 1 2 ^(0^ — 2) 1 ^ \ 

^J^(^2 + M2 + ^2)a+2 - r(a + 2) (/.2)a-2i ^ y-2(a-2)i^j | 

(A12) 

with 

r ^ r ^ N\ r 



where 

Mf = /i^ + /M^ , = /i„ (AT) In N (A14) 

^{fiupa} = S^uSpa + S^pSua + ^ pa^ pu (Al6) 

( , y ^ n( x/{lK{N)\nN) 
l/-.(.-.)(x) = -limX.q (^ i^,/(^;,^(^)i,^) J (Air) 

7^ = 7^ = 0.5772 ■ ■ ■ , /i^(iV ^ oo) = 1 (A18) 

The exphcit form of y-2{a-2){.x) is aheady given by Eq.(l9]) 

Appendix B: Details of Calculation of Self-Energy Insertion Diagram 

In this Appendix, we shall present the details of the calculation of the UVDP parameter 
integrals for the self-energy correction insertion diagram (ai) shown in Fig. ([2]). 

With the simphfication by setting /i^ ^ and yi{-^) — )• in Eq.(l2T]). we have the 

following expression for Ai 



K) 

2 



' ~ (levrT Jo 1 + Wo Jo il + v^r6' 

= Mil^^ + M^^'^ + Mi^^l (Bl) 

Here A^22^'', -^20 ^ and M.'^^^ represent quadratical divergent, logarithmical divergent, and 
regular part in the Ai^2^\ 



^ - '-^-rMi I'd. f-^ 

22 (167r2)2 6^ l + uJo Jo {l + vi^ 

(167r2)2^ ^2 Jo 1 + u 



'Jo ql-ql 

where in the third line we use the trick that = —q'^u, which can effectively transform a 
UVDP parameter integral into a 1-folded irreducible loop integral (ILI) that is the funda- 
mental object to be regulated in the framework of LORE. In the following, we will frequently 



use this trick without any explanation. 



-P^ 2 Jo 



(167r2)2" " ^ -p2 ' ^2 7o (l + t^i)=' 2 3(l + t;i) 



(167r2)2^ ^ ' -p2 ' ^ ^6 7o 1 + ^1 36 



(167r2)2^ ^ ^ -p2 ' ^ ^6 7o g?-go 36' 

(^^y (In ^ _ + 1) . [ (In ^ - 7,) - — ] (B3) 



(167r2)2^ ^ ' -p2 ' / IgV 36 

The integral 



_ 8.e^r(4) l rfnn 
■^^^ - (16vr2)26^ ^ ^Vo {l + ufjr j. il + v,r 

(B4) 

is difficult since the complicated form in the logarithm. However, if we only focus on its 
divergence behavior, we can make a lot of simplification. As discussed in Section 2 in 
our previous paper [2|, the analogy between circuit diagram and Fig. ([2]) tells us that the 
divergence appear only in the region vi — )■ oo. Therefore, we can set Vi in the region vi > V, 
which \^ is a very large number V ^ 1. In this asymptotic region, we have — t- 1 and 
the above expression simplifies to: 

dvi 



[In ^ + In (x{l - x)] + ln(l + Vi)] 

[1 + uY V / 
^ Mt^^ + Mt^ + M^i^, (B5) 

where M.^2Ru^ -^2iji ^^"^ -^2'r»!i represent the three parts in the second line, respectively, 
the results of which are given as follows: 

Mtl = ■7^^{-a>"p') / '"TT-^ / dxx{l-x 



sfe" , „„ ,„ ,,1 r A'. 



j'V)(-i)iX r 



(167r2)2^ ^ ^ 'Qj^lJ^y^ 

^(^7^y)(ln , ^ -7j, (B6) 



(16772)2 3'" " -g2(l + y) 



1 + t;i 



8^e^ . -^-V)i(-A) 



(167r2)2^ ^ ^ ^2' 18' Jv- 1 + ^1 

VV(ln^#^-7.), (B7) 



(167r2)2l8" ^ ' -ql{l + V) 



ja,) 8ie*^ , o, duu /"^ , r M /"°° dvi 



Since we are only interested in the divergence behavior, the last finite term can be neglected 
due to the complexity of its calculation. Below we only give the regulated result for the first 
term with LORE method and leave the details of derivation to Appendix D. 

-^?'^V)[(ln^-7.)^-«.]. (B9) 



Putting Eq. ( ]B6[) . ( 1B7I) and ( 1B9I) together, we arrive at the result for Ai'^^^ 



Mir ~ (IS^S-'Vl^dn - 7.) - ^(in^l - -,.y + -aj. (BIO) 

With the experience of the calculation of Ai^2^\ the computation of logarithmic part 
A^o"^'* is straightforward: 



(167r2)2 6 Jo {l + uyjQ Jo {l + viy l + vi 



M2 a;t;i 
^ A^i;^)+Al("^^\ (Bll) 



^From the definition of A^oo(/?)' ^^^^ ^'^^^ ^'^^y essentially the same as 

counterparts in the quadratic part A^20(i?)- Thus, we can write down the final expression 
for A^oo(i?) imniediately following A^2o{i?)' 

By summing Eqs. (IB2p . (]B3p . (IBlOp . (lB12p . (]B13p . we finally obtain our result for diagram 
(ai) 

(167r^)2 — 6 — g2 

1 M2 A/f2 1 m2 1 

+ (^My _ pMp.) . [ (In ^ _ ^ _ _ (In ^ _ ^J2 ^ 

5 , 11,, 



Appendix C: Details of Calculation of Vertex Correction Insertion Diagram 



This appendix gives the details of the calculation of integration of UVDP parameters for 
the vertex correction insertion diagram (6). 

Note that the integration of k2 naturally separate M.^''^ into four pieces with different 
powers of factor ^ i — , each of which we denote as M.f\ i — 0,1,2,3. In each 

""'"(l+«'l)(l+"2) 

piece, there are different degree of divergences, quadratically or logarithmically divergent 
and regular. In the following, we will use a second subscript 2, 0, to represent these parts. 
Let us compute these parts one by one. 

Firstly, M.q^ contains only the single quadratically divergent term M.q2-> whose result 
can be given directly as: 

(,) _ 8^e^m r-rr T ' ).-'j^(0) 

- (167r2)2 36^ ^^^> 

where, as noted before, we have set — )■ in the end since it can be regarded as just the 
IR regularization in the LORE and does not affect the UV structure of Feynman integrals. 
In the following, we will always apply this method to simplify our results. 

Also, M-f^ merely includes one term M.12, but the calculation is a little involved since it 
can be divided further to several pieces A4i22, •^120 •^i2ii according to the degree of 
divergences each of them contain: 

,1 2 



u + 



(l+Vl){l+V2) 



_ 8ie^r(4) rn t ' ) r - 

1\ Jo l\{i+v^r t^i+^/io 



{^c-/^^[ln^-7. + l]} 

= Mfl + Mfl + Mfi^. (C2) 
Each term above can be calculated straightforwardly: 



[ibn ) Jo Jo i^^[i- + Vi) i^^i- + v^Jo ii+ 



i=l ^ 1=1 



(l+t,i)(l+„2) 

-^eb(ln-|-7a.) + T^], (C3) 



(167r2)2 ^^6' -g2 IS' 



1 r 1 r (a;2-a;i)^ ^ a:i(l - xi) ^ X2{1 - X2)^ 



(l + ^i)(l + ^2) 1 + ^1 1 + ^2 

1 (X2 - 



[^+ (1H-..)U.) ]M1 + ^1)^(1 + ^2)- 



8ie^r(4) 2 13. M2 

V^(ln— |-7u,)- (C4) 



Since it is easy to prove that Aif2R finite and thus irrelevant to our present discussion of 
UV divergence structure, we omit its calculation here and simply write down the result for 
Mg as: 

8ie4r(4) 1 , ^ 5,. 13 2n 

where ~ means that they are equal up to the divergence part. 

By observation, Aif^ contains quadratically and logarithmically divergent parts, A^g 
and A^2o- The computation of A^g is much like that of A^g and the result is shown 
below: 

= Alfi + Mgo + -Mgi. (C6) 

8^e4r(4) /-^A. r A ... ^ 1 



_ _8^e^ r-p. r-TT dv, 
•^222 - (ig^2)2 ^ X fj Jo l\ (1 + -.)^^ h 1 



00 

du 



1 



(l+l>l)(l+D2) _^2 



N + (l+^,i)(l+^2)]^ 



_8ie4r(4) 

220 - ~ nR-^2N2 



io l\ Vo fi(l + ^.)^ + Jo 



(i+«iKi+f2) r r (^2-a;i)^ ^ - xi) ^ 3:2(1 - ^2)^ 



[«+(TWaT^P(i + -i)^(i + -2)^^^ -p' 



(167r2)2^ ^ 36-6^ -p2 ^'"^^ 



We can prove that the term A^22i? is finite. So we have the divergence behavior of the 
integration is: 

AI20 contains many different Lorentz structures, so it is useful to further refine it into 
three parts according to their Lorentz structure: represents terms proportional to 

p'^p'^, M202 proportional to g'^'^p'^ and M^os proportional to g*^^ 

_ 8ze^r(4) r^^ fn 5(1 V ^ ) Hdu 

= A^Sio + ^Siii, (CIO) 
where the last subscript of Al2oio(i?) represents the divergence degree of the part. 



8ie''r(4)17 „ 

As before, we can prove that the part Al2oii? is finite, so the divergent part for M.201 is: 

By the very similar way, we can obtain the result of A4fo2- 
1 ri 6 1 {x2-xif 

(l+fl)(l+^'2)J(l+^l)2(l+t;2)2 



+ :;- o)(- 



l + Vi I + V2 2''1 + Vi I + V2 {l+Vi){l + V2)u + 



(l+i;i)(l+«;2) 



2 



^1 , ^2 w, Xi X2 



l + Vi I + V2 l + Vi I + V2 11 



^^0^20 + M)-^^ (C13) 



By integrating the expression before j^{ln — 7^^), we obtain M.2Q2q 

The part of A^202K can also be proved to be finite, so the divergence part of M.2Q2 

^Jb) 8?e^r(4) 5 „^ 2/1 -^c N /r.-.rX 

Much like X^, and Af^i 

in structure, it is expected that the A4203 ^e calculated 
similarly. However, a very important new feature will appear: the overlapping divergence 
structure hidden in our current LORE procedure. 

(,) _ SieMl) 2 f A^^. rn J(l V M Tc^u 

•^203 - (ig^2)2 5 ^yofi ^ h^ + '^^ Jo 

1 . 1 1 , 

+ (TT^MTTT^]' 6+(l + .0(l + -2)^' 
(■ 1 r (x2-a;i)^ _^ - Xi) _^ ^2(1 -^2) . 

^+(TT^;Ty(TT^ (1 + ^0(1 + ^2) l + I + V2 

1 (X2 - 



[^+(TT^;rfei)]Mi + -iRi + -2> 

([ln-^-7.]- 



^^1- 1 r (3^2-3:1)^ _^ Xl{l - Xl) _^ X2(l - X2) . 



1 {X2 - Xl) 



[^+(iT^;iiT^]Mi + -im + -2)^ 
= ML + ML, (CI6) 

where, A42030 represents the part proportional to (In ^ — 7;^) while A4203i? ^^^^ 



regular terms. The integration of A^2030 straightforward. 



(,) _ 8.e^r(4) , r^f. fn 5(1 T ' ) Hdu 



1 1 1 M. 



2 



1 ^ (a;2-a;i)2 , a;i(l - a^i) _^ X2(l - ^2)1 



1 (a;2-a;i)^ , 



— — ^c/^ p In — 2~^^nT^-^y — +/ ^ — 

(167r^)^ -p^ 12 72 ^0 1 + f 1 Jo 1 + ^^2 

W(ln— |-7-)[T7^-7T7T(ln-|-7-)]- (C17) 



(167r2)2 ^ ^ ^ -p2 '-^^12 36^ -g2 

Note that the above parameter integration before (In — 7^;) is also divergent, which 
is the signal of overlapping divergences. It is clearer when we use the analogy between 
the Feynman diagrams and the electric circuits, which tells us that the integration of Vi 
and V2 should reproduce the divergences coming from the subdiagrams of left and right 
vertex corrections in Fig.dl]) since (1+^) times the effective propagator kf + xi{l — 
Xi)p'^{[{ki + k2 + {x2—xi)p)]'^ + X2{l — X2)p'^) for the left (right) circle in Eg. fl26|) . Whent>i(t>2) 
tends to infinity, the left effective propagator ^^^[^1 + ~ Xi)p'^] (the right counterpart 
j^iii^i + ^2 + {^2 — Xi)p)Y + 3^2(1 — 2:2)^^}) approaches zero, which means that the 
left (right) half circle collapse to a point and that sub-circuit is short-cut in the electric 
circuit language. This singular behavior will become manifest as the divergence in the final 
integration of UVDP parameter fi(f2), just as the ones shown in Eq.( IC17p above. Thus, 



the result in Eq. flC17p can be understood as the sub divergences coming from left and right 
vertex correction times the overall one (In — 7^^), which is the definition of overlapping 
divergence jl, 0] ■ 

This overlapping divergence take us some further difficulties that the integration of A^203R 
turn out to be divergent, which gives us further contributions to our UV- divergence structure. 



1 r I 1 

+ 



I 1 (a^2-a^i)^ _^ a;i(l - xi) _^ 3:2(1 - X2) . 

^^+ (l+,,)(l+^,) Hl+t^l)(l+t^2) 1 + ^1 1 + 1^2 
1 (X2-Xi)^ 

^ r (3^2 ^ - xi) ^ 3:2(1 - X2) 

^+ (i+.oW^.) + + 1 + ^1 1 + ^2 

1 (3^2-3:1)^ 



r}- (C18) 



Because of the complication in the logarithmic function, it is difficult to get a closed an- 
alytical expression for this kinds of integration. However, if we only focus the divergence 
behavior of the integral, then we can use the method introduced in the calculation of M.'^^^ 
to greatly simplify the integral and to obtain the asymptotical results. To disentangle the 
overlapping divergences, we need first to know in what parameter space region the diver- 
gence happen. Prom our experience when working M-foso, the divergences take place when 
Vi — )■ 00 and V2 — )■ 00. 

Now we consider the region where vi ^ 00 and V2 ^ 0. Fist we choose a very large 
number, say V, and set the integration region is only confined in V. In such a region, 
are small quantities, so we can expand the expression according to The leading 

term of M.203R is. 



^ 8.e>r(4) 1 r^f. r- dv, r- ^^^^{i^ X2{i-X2) 

■ Ig'^p'i f dX2X2{l - X2) HX2{1 X2)] T H du \ 

(167r2)2 6 Jo Jv (l + 'yi)^7o {u+t^Y 



- [ (/a;2X2(l -X2) / [ du- — \ ln[M + 

Jo Jv i^ + vi)^Jo («+iT^)^ l + ^^i 



-]} 



8^6^1(4) 1 2r/ r dvi I A r dvi 1 p dvi , 

r9^P{{—rz)^ ^ ^ 0/ ln(l + ^^i)]} 



8ie^r(4) 1 „^ 2rl3„ M2 ^ 1 , 1 . .^^ , 



(167r2)2 72" "^6^ -ql{l + V) 2' -ql '"^ 2 

where in the last line we extended the lower bound of the integration range to as before 

for convenience. 

Since in the other asymptotic region t;2 — )■ 00, — )■ 0, we can obtain a similar expression 
except for the exchange of vi •<-> V2 and 3:2 •<-> xi. Thus, we can expect to get the same 



asymptotic result. Therefore, the divergence behavior of A^203i? 

8ie^r(4) 1 „^ 2rl3,, Ml , 1,, ,3 1 ^r^nnx 

Alg^^ can be naturally divided into three parts AI32', AI30 and Alg^ according to the 
divergence degree followed by integration of loop momentum k2- M!"^2 contains only one 
term, so the calculation is straightforward: 

_ 8ze^r(4) 1 r^^ fn 5il T ^ ) Hdu 



(l+i;i)(l+t)2) 



^^2^2 + ^^2^0 + M^^^ (C21) 



1 



[«+(TT^;Tfc^F(l + -im + -2)^ 

(167r2)2 36^ ' ^ ' 



1 1 , 

a + ^i)(l + ^2) l + l + ^;2 

1 (3:2 - Xl)2 

[^+(iT^;Tfc^]Mi + -im + -2)^^ 

bmce we can prove that Mf^^ is finite, the divergence structure we concern now is given by 
the addition of A^322 and AI320: 

Like M^2o-> -^30 contains three parts differentiating by their Lorentz structures: Aifoi rep- 
resents part proportional to p^p", Mf^i for g'^'^p'^ and AI303 for g^^/J^^- Let us first calculate 



MS- 



1 
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= M^L + AlSiR- (C25) 

Actual calculation of A^3oio shows that the integral before (In ^ — 7^) is logarithmically 
divergent. So like the part A^2030; -^3010 ^^^^ involves overlapping divergences. 
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For the more challenging part M-fQi^, we follow the approach already used in the deriva- 
tion of A^203K order to obtain only the asymptotic results. We first focus on the region 
f 1 — !■ 00 while t>2 — )■ 0, the expression for A^3oi/j can be simplified to: 

(C27) 

where the integration over vi is confined in the region from ^ ^ 1 to infinity. It is easy to 
see that the simplified integral is essentially the same as that for A^203ii' (EH]), so we can 
straightforwardly write down the result for A^3oir 
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where we have already doubled the result for A^3oi_Rt,i due to another asymptotic limit region 
V2 — 7- oo,t>i — )■ 0, which gives the exactly same result. 



A very similar calculation as Aifgi can give us the results for M.fQ2, 
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Note that the coefficient before (In ^ — 7^) is divergent, which means that M.^^q2 ^^^o 



contains overlapping divergences. In order to deal with the resultant divergent integral 
^to2Ri following the procedure used 
integral in the region i>i — >■ 00, 1)2 — >■ 



^to2Ri following the procedure used by M.2mR ^^"^ ^toiR^ ^® need first to simplify the 
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In order to obtain A^302R! the above result needs to be doubled to take into account another 
contribution from the region V2 00, Vi 0, which is given below: 
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For the accomplishment of the computation of M^q, we have to calculate Aif^^: 
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= Al(?3o + A^Sr, (C33) 

where A^soso represents the part proportional to the logarithmic divergence coming from 
while A^303/j the rest parts. It is a direct exercise to obtain A^303o and the result is: 
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while part A^sosr can be proven finite. So the divergence structure of Aif^^ is given by: 
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For the part M.f]i, we can prove that most terms are finite for the integration of UVDP 
or Feynman parameters, except for the following one: 
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By putting the divergence parts of the terms Eqs.dCl]), flCT2|) . f lCTS]) . flCTTl) . 

( 1U20D . ( Oil) . ( Oej) . ( OSj) . f OOj) . ( 02]) . f lUSSj) . (ESg) together, we finally arrive at the 



divergence behavior of the diagram 
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Appendix D: Derivation of the Integration Jq°° in(l + V2) 



Recall that in the above derivation of the two-loop massless QED photon vacuum polar- 
ization diagrams we encountered a new parameter integral 
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which does not appear at one-loop level. Thus, for completion of our calculation, in this 
section we shall derive its regulated result with the LORE method. 

With the same philosophy as before, we would like to transform this UVDP parameter 
integral to a momentum-like one. In order to do so, we just need to multiply a free energy 
scale — in the numerator and denominator simultaneously, 
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= 7i - 72, (D2) 

where we have defined q^ = —q^V2 and separated the integral into two parts /i and 12- The 
integration I2 can be easily worked out with the LORE method since it is exactly the one 
encountered in the one-loop calculations. The result of I2 is 

72 = ln(-0[ln-f-7. + yo(^)] (D3) 

where wc define fJ'l = fi'i — qo- 

The integration 7i is really new and requires us to calculate more carefully. In the 
following, we shall give the detailed derivation of this integration in the LORE method. As 
usually done in the LORE method, we need to apply a series of regulators to the integral: 
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where we defined = —ql + jj^^ + lM'^ = jj^^ + lM'^. With the help of the following equality: 
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we can easily work out the above integration: 
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where we have introduced two series: 
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both of which in the hmit of — ?■ oo can be shown [3] to have finite hmits 

Ml = hm Ml{N), and 7^ = hm 7^ (AT). (D8) 

We also define a new function 
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since the expansion of ln(l + is of the order of 0{-^). The constant a^j is defined as 



which can be seen that when the UV scale tends to infinity the function y'o{j^) vanishes 
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It can be shown from our numerical calculation that a^^ is finite and is = 1.62931.... 
By combining the expression Ji and I2, we can give our final regularized result 
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= 2 [(In -f - l^f + 2 In ^(In - 7.) - + 2 In fJl/o(^) + l/;(i^)].(Dll) 

When we set the IR scale /i^ to which can be viewed as the IR cutoff in the LORE method, 
then the regulated integral can be simplifies to: 



If we further take the hmit Mc — )■ oo, the last term yoi^r) — )• also, which gives the result 
we use in our previous calculations. 
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